Intersubband Polaritons in the Electrical Dipole Gauge 
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We provide a theoretical description for the coupling between the intersubband excitations of a 
bi-dimensional electron gas with the electromagnetic field. This description, based on the electrical 
dipole gauge, applies to an arbitrary quantum heterostructure embedded in a general multilayered 
waveguide or a microcavity. We show that the dipole gauge Hamiltonian automatically takes into 
account the Coulomb interactions in this system. Furthermore, it can be conveniently expressed in 
terms of the many-body collective plasmon modes, which interact both with each other and with 
the light field. The dipole gauge therefore provides a suitable framework for the study of solid state 
Quantum Electrodynamics (QED) phenomena, such as the ultra-strong light-matter interaction 
regime, occurring at very high electronic densities. 
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I. INTRODUCTION 

In the description of a physical problem dealing with 
light-matter interaction there is always a degree of free- 
dom in the choice of the potentials associated with the 
electromagnetic field. Even though the physical phenom- 
ena are obviously independent from the particular vector 
or scalar potentials used, the choice of the gauge is a 
critical issue, as it can be more or less adequate for the 
physical interpretation of the phenomena. In the liter- 
ature there is already a consensus that the interaction 
of non-relativistic bound charges is very conveniently ex- 
pressed in terms of the dipole gauge This has been 
firstly identified by Power and Zienau in the 50 's 
and subsequently by Woolley In this formulation, 
called Power-Zienau- Woolley (PZW) gauge, the sources 
are electric and magnetic polarization fields and the cou- 
pling occurs via the intensities of the displacement field 
D and the magnetic field H rather than using the scalar 
and vector potentials V and A. In the case where the 
magnetic interactions in the system can be neglected, the 
PZW gauge is also called dipole gauge since the interac- 
tion Hamiltonian contains only the coupling between the 
material polarization P with the light field. 

In our article we propose to apply the dipole gauge to 
the case of intersubband transitions interacting with a 
photonic cavity mode. The motivation of our approach 
stems from the fact that we are interested to investigate 
a regime, called ultra-strong coupling firstly introduced 
by Ciuti et al. [1]. This regime of light-matter interac- 
tion is attainable in a bi-dimensional electron gas with 
very high density embedded into a photonic microcav- 
ity and is characterized by the fact that the cou- 
pling and the material excitation energies are comparable 
quantities. However, the very high electronic densities in- 
crease on one hand the light-matter interaction and on 
the other renormalise the transition energies of the sys- 
tem due to a collective effect, the "depolarization shift" 
Q • This effect cannot be neglected in the limit of ultra- 
strong coupling, however it was not explicitly considered 



in the initial study described in Ref. 4]. In our work 
we show that the dipole gauge Hamiltonian handles the 
light-matter interaction and the depolarization effect on 
the same footing. The underlying physical picture is that 
the active polarization that couples to the photonic mode 
is not a single particle electronic transition between con- 
fined states but rather a collective electronic mode, - a 
plasmonic mode -, arising from electrons distributed in 
different subbands, yet phased by the Coulomb interac- 
tion. The latter point is one of the main conclusions of 
our theoretical investigation. 

The correct description of the collective electronic ex- 
citations and their interaction with light can be obtained 
in the Coulomb gauge, only if both the vector potential 
A and scalar potential V are considered. This has al- 
ready been noticed in studies of ensembles of two level 
systems interacting with light 0. In the absence of the 
cavity, thus in the weak coupling limit, collective effects 
can be identified with the scalar potential V @, which 
describes an instantaneous interaction between electrons 
p^ . However, when dealing with the resonant coupling 
of a microcavity with an electronic transition the ob- 
servable quantity is the retarded electromagnetic field 
that can be obtained through a proper combination of 
both V and A 10]. In order to have a general Hamil- 
tonian, which is valid from the weak to the ultra-strong 
coupling regimes, it is therefore essential to include the 
Coulomb potential V, thus completing the study of Ref. 

Such Hamiltonian is readily obtained in the dipole 
gauge, where the matter degrees of freedom are gathered 
in the polarization field P, which not only describes the 
interactions between the electrons, but also couples to 
a material independent and retarded photonic field de- 
scribed by the electric displacement D [l|, . Indeed, as 
we previously pointed out in Ref. @ , the depolarization 
effect in the bi-dimensional electron gas is contained in 
the quadratic term of the dipolar Hamiltonian. As a 
result, our model correctly describes the local field effects 
arising from the very different spatial scales between the 
electronic and photonic confinement. The results of Ref. 
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|j| are contained in our formalism and recovered only in 
the opposite limit, i.e. when the electronic polarization 
fills the whole cavity volume. The importance of the spa- 
tial overlap factor was also outlined in our experimental 
study of the ultra-strong coupling regime [y, as it 
allows to correlate both polaritonic and weak-coupling 
absorption data. On a more fundamental level, we show 
how the spatial confinement of the electronic polarization 
leads to the No-go theorem for intersubband transitions 

Our paper is organized as follows. In part|TT]we estab- 
lish the Hamiltonian of the system in the dipole gauge. 
We consider a very general case of an arbitrary het- 
erostructure embedded into a general planar waveguide 
multilayer. The microscopic expression of the intersub- 
band polarization field is derived in section IIIDl in the 
long-wavelength approximation. This microscopic ex- 
pression relies directly on the electronic wavefunctions, 
and therefore allows to go beyond the semi-classical 
model employed in Ref.Q, as it applies to an arbitrary 
heterostructure potential, with an arbitrary number of 
occupied subbands. 

In part IIIII we adapt the general dipolar Hamiltonian 
obtained in part|lT]to the case of a series of highly doped 
quantum wells and we express it in terms of the collective 
electronic excitations of the system. We also provide a 
version of this " plasma Hamiltonian" for the case of OD 
resonators fsection IlII CI) . In section UlI Dl wc study the 
correspondence between the plasma Hamiltonian, trun- 
cated for the case of a single intersubband transition and 
a single waveguide mode, and the corresponding Hamil- 
tonian in the Coulomb gauge. This allows us to connect 
our formalism with previous work In particular, we 
show that the dipole gauge provides automatically the 
relevant contributions due to the Coulomb interaction of 
the system, that were missing in Ref. [ij . 

The formalism is applied to study the properties of 
the polariton states in part IIVI In section IIV Al we ex- 
amine the polariton dispersion, and in section IIVBI we 
discuss the No-go theorem for intersubband transitions. 
We show in section HV CI that our formalism is consistent 
with the effective medium approach, in full agreement 
with what has been proposed in the literature IJ] . Most 
of the technical details have been gathered in the ap- 
pendixes. 



II. INTERACTION HAMILTONIAN 

A. General considerations 

The interaction between light and quantum het- 
erostructures is usually studied in planar multilayered 
systems. This geometry is most naturally compatible 
with the epitaxial growth. A very general system is de- 
scribed in Figure [ija). It consists of homogeneous, non- 
absorbing and non-dispersive dielectric layers described 
by real dielectric constants e^, embedded between two 



infinite semiplanes that act as optical claddings. We 
can then define a piecewise dielectric function e{z), with 
e{z) — £i in the layer, z being the growth axis. The 
multi-layered dielectric stack defines guided modes that 
confine the light field around the heterostructures. Note 
that, at this point, e{z) does not include the resonant 
contribution from the electronic transition in the het- 
erostructres. This contribution will be included through 
the coupling between the guided modes and the electronic 
polarization. The full Hamiltonian of the system writes: 



H — Hp + Hn 



(1) 



Here He is the Hamiltonian of the electron gas in the 
heterostructures, which is provided explicitly in section 
III CI and Hp is the photon Hamiltonian: 
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with D(r) and H(r) respectively the displacement field 
and the magnetic field operators. The interaction Hamil- 
tonian _ffint in the electrical dipole gaug e is written as, 
neglecting the magnetic interactions [15| : 



H, 
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eoe{z) 



D(r).P(r) + ip2(r)' 



(3) 



Here P(r) is the polarization density operator of the 
electron gas. This is a central quantity in our theory, 
and it is provided explicitly in section IIIDl from a micro- 
scopic model. One issue that we will discuss in details in 
the article is the role of the quadratic interaction term, 
P^(r), which describes the self- interaction of the elec- 
tronic polarization, and therefore contains the effects of 
the dipole-dipole interactions. 

Most generally, the expression of the interaction (|3]) 
should be written as a non-local expression, which takes 
into account the spatial dispersion of the electromagnetic 
response of the medium [16j. However, as illustrated 
in Figure (Ua), in the systems that we study, the het- 
erostructures embedded in the multilayered stack have 
a typical extension in the growth axis (the z-axis) that 
is much smaller than the wavelength. This justifies the 
local form of the interaction postulated in ^ . 

We can evaluate the convenience of the dipole gauge 
already at the classical level, where the displacement field 
D(r) is determined only by the charges exterior to the 
system. On a quantum mechanical level, this means that 
D(r) describes a purely transverse field Therefore, 
neglecting the dissipation in the system, we can consider 
D(r) as a free photon field, independent from the elec- 
tronic polarization of the heterostructures. This can be 
clearly observed, for instance, in the normal component 
of the displacement field Dz, that, while it couples to the 
intersubband polarization, is continuous across the inter- 
faces. We therefore consider the derivation of D(r) as a 
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FIG. 1: a) General planar multi-layered system with a piece- 
wise dielectric function e{z), supporting guided modes. There 
are quantum heterostructures embedded inside the multi- 
layered stack, b) TMo mode guided between two metallic 
plates. 



separate problem, which is dependent only on the partic- 
ular arrangement of the multilayered stack. In this ap- 
proach the resonant contribution of the electronic inter- 
subband transitions taking place in the active media are 
contained separately in the polarization operator density 
P, and are active only through the interaction term iJint. 
Later on (section lIVCp we shall see that this approach 
leads to an effective medium treatment of the system. 

In our treatment the dissipation effects will be ne- 
glected. The latter can be taken into account as in the 
classical paper by Huttner and Barnet [l7| . by including 
a dissipative bath with continuous degrees of freedom. 
This approach will however require more precise defini- 
tion of the displacement field D, as it leads to a noise 
contribution from the dissipative bath [l^, • Here we 
will restrict to a fully Hamiltonian treatment. The dissi- 
pation will be taken into account only in a phenomeno- 
logical way in section IIV C[ by adding a small imaginary 
part to the eigenfrequencies of the system [Tsj . 



B. Free photon Hamiltonian 

For the quantum description of electromagnetic field, 
we use the basis of the guided modes of the multilayered 
stack, that are bounded in space. Due to the transla- 
tional invariance of the system in the plane perpendic- 
ular to the z-axis, the guided modes are characterized 
by their in-plane wavevector q and their energy fiwcq- 
The function — Wcqdql) defines the dispersion rela- 
tion of the guided modes. The heterostructures interact 
only with the TM-polarized modes, in order to respect 
the selection rule of intersubband transitions [l^. Us- 
ing the general quantization rules HO] , we can assign 
bosonic creation and annihilation operators a]^ and Oq 
to each guided mode. The components of the quantized 
electromagnetic TM free field are then: 



H = ^ zAq(e, A e,)/q(z)e''i'-|| ^ „t_^) (4) 
q 

q 

D|| = - ^ Aqeq-^^e-ii (a, - al^) (6) 




Here S is the area of the system, and r|| is the in-plane 
position vector. We have introduced the unit vectors 
and Gq = q/|q|, and the symbol "A" designs the vector 
product. The constant Aq is the vacuum field intensity of 
the guided modes. It has been derived, for this particular 
system, in Appendix[Al The function fq{z) describes the 
lateral profile of the guided modes, and Lq is a normal- 
ization coefficient: 

/+00 
Uzfdz = Lq (8) 
-00 

If we set the maximum of the dimcnsionless function 
/q(z) equal to 1, with this definition Lq can be regarded 
as the effective thickness of the multilayer system. The 
function fq{z) satisfies the Helmholtz equation: 

According to the boundary conditions for the electro- 
magnetic field, /q and £{z)~^dfq/dz must be continuous 
at the interfaces between the different layers, and fqiz) 
must vanish for z — > ±00. Such a function has been illus- 
trated in Figure [IJa). Equation © will have in general 
multiple solutions, describing different modes labelled by 
a discrete index j, that can be added in the notations if 
necessary. 

One or two boundaries of the multilayer system can be 
metallic. In the mid- and far- infrared frequency range, 
the metals are described by a very large negative real 
dielectric constant em < 0, and the electromagnetic field 
density in the metallic layers is vanishing. Therefore, in 
order to have a consistent quantification scheme for the 
electromagnetic field, without the burden of quantifying 
the conducting electrons in the metallic boundaries, we 
neglect the electromagnetic density energy in the metallic 
regions in the integrals ([2]) and ([3]). However, the finite 
dielectric constant sm can be taken into account through 
the dispersion relation of the guided modes, as for the 
example provided in Appendix [Xj 

With the use of the above expressions, the Hamiltonian 
([2]) of the free electromagnetic field takes the standard 
form: 
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FIG. 2: (a) Typical quantum heterostructure potential V{z): 
a quantum well. There a two bound levels with wavefunc- 
tions (/>i,2(z) (b) General potential of a doped heterostruc- 
ture, including the Hartree correction Vh(z) due to the static 
Coulomb interactions. 
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(10) 



We now apply this formalism to a particular system. 
For mid-IR and THz frequencies the simplest confining 
system is provided by the TMq mode guided between 
two metallic plates, separated by a distance Lcav, as il- 
lustrated in Fig. [ijb). The semiconductor with a di- 
electric constant e fills the space between the plates. In 
this case, in the limit of a perfect metallic boundaries we 
have fq{z) = 1 and Lq = Lcav and the dispersion rela- 
tion becomes = ew^^/c^ (Fig. mb)). The in-plane 

component of the displacement field vanishes, Dy = 
and the remaining non-zero z-component is: 



D, = ie. 



E 

q 



(11) 



In the rest of the paper, we consider exclusively the 
TMq mode, which simplifies greatly the calculations 
without a loss of generality. In section IIII CI we also con- 
sider the case of double metal zero-dimensional (OD) mi- 
crocavities in which we have added a lateral confinement. 
In these structures the propagation wavevector q of the 
TMq becomes quantized. 



C. Electronic Hamiltonian 

In a semiconductor heterostructure the band offsets 
between different semiconductor layers provide a confin- 
ing potential V{z) in the growth axis. Typical exam- 
ple of a confining heterostructure, the quantum well, is 



sketched in Figure ^a) , and a general heterostructure is 
illustrated in Figure [H^b). We consider the case where 
the Fermi level Ep lies into the conduction band, due 
to intentional doping, which leads to the formation of a 
bi-dimensional electron gas [7] (figure [Sfb)). The poten- 
tial V{z) yields discrete energy levels fujj\, labelled by 
an integer index A. Electrons are free to move in the 
plane perpendicular to the growth axis, and we denote 
by hh their in-plane momentum. This free movement is 
described by the parabolic subbands illustrated in Fig- 
ures [2Ja) andm^b), and the total energy of an electron in 
the subband A is: 



2m* 



(12) 



with m* the effective electron mas. 

In a highly doped heterostructure, in order to deter- 
mine the confining energies fuvx, one must take into ac- 
count not only the heterostructure potential, but also 
the Coulomb interaction between the charges. In the 
Hartree approximation, this static Coulomb interactions 
are described by a self-consistent potential Vh{z) due to 
the presence of electrons and ionised impurities 0, fioj . 
The potential Vh{z) depends on the envelope wave- 
functions (j)x{z). Therefore the energies hux and wave- 
function ipxiz) are determined altogether by solving the 
one-particle Schrodinger equation with a total potential 
V(z) + Vh{z) self-consistently with a Poisson problem Q: 



2m* dz^ 



V{z) + VH{z)\Mz)^f^xM^^) (13) 
d^Vniz) 



dz^ 



P(^) 



e 

££o 
m 



Nd(z) 



(14) 



^7V.|0.(z)P (15) 



Here p{z) the electronic density, Nd{z) is the dopant 
density, and Nx is the population of the A-th subband. 
The exchange-correlation effect has been neglected in the 
above equations. This set of equations leads to one- 
particle quantum states |A, k) with normalized wavefunc- 
tions: 



(r|A,k) = 0A(2)-^exp(ikr|| 



(16) 



The corresponding fermionic destruction and creation 
operators are c^k and c^j^. The electronic Hamiltonian 
acquires the one-particle form: 



i7e = ^ ftwAkC^kCAk 



(17) 



Ak 



At this point we have determined the stationary state 
of the system. The effects of the Coulomb interaction 
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that we have considered so far are the static effects 
arising from the inhomogeneous spatial distribution of 
charges. These static effects have been lumped into the 
one-particle subband energies hioxi^, through the self- 
consistent set of equation (13)-(15). 

In the next section we use the stationary basis of single 
particle states (/'a(-z) in order to define the microscopic 
polarization density of the electron gas. This will allow 
us to study the excited collective states, coupled with 
light, and to recover the dynamic effects of the Coulomb 
interaction, such as the depolarization shift. 



D. Microscopic expression of the Polarization 

Usually, for the studies of the interaction of intersub- 
band transitions with light, one does not define a polar- 
ization operator as an independent local quantity P(r), 
but rather computes the linear response of the electronic 
system due to the solicitation of an external harmonic 
electric field. The linear response is described by the fre- 
quency dependent non-local susceptibility x('^)i211] which 
is computed from the current-current correlation function 
through the Kubo formula ^3 ■ The collective excitations 
of the system are then obtained from the isolated poles 
of xi^)- 

In order to use the electrical dipole gauge, as formu- 
lated by the interaction Hamiltonian Q, we need to ex- 
press the local polarization operator P(r) as an indepen- 
dent quantity. Classically, the local polarization P(r) is 
defined as the average dipole moment of the charge distri- 
bution over some microscopic volume (23j . For instance, 
in the case of the square quantum well with a thickness 
Lqw illustrated in Figure HJa), such volume could be the 
volume of the quantum well. However, it is difficult to 
apply such a definition for an arbitrary heterostructure, 
such as the one depicted in Figure [IJb), where the spa- 
tial extension of the confinement potential varies with 
the energy. In this case, the averaging volume becomes 
an arbitrary quantity. Since the microscopic intersub- 
band dipole is defined from the electronic wavefunctions 
(j)\{z), the truly microscopic expression of P(r) should 
involve directly 0a (^). 

Similar problem has been encountered in the attempts 
to define the static polarization in ferromagnetic materi- 
als as function of the electronic wavefunctions and nuclei 
distributions [13]. Then the idea is to define the polar- 
ization as the time integral of the microscopic current 
arising during the adiabatic switch from one configura- 
tion of electrons and nuclei to another. 

In our case, the microscopic current which corresponds 
to the intersubband transitions is a rapidly oscillating 
function at the frequencies of the transitions. Then we 
define the polarization in such a way that its time evo- 
lution under the full Hamiltonian would lead to a micro- 
scopic current: 



^ = -[P(r),7?]=j(r) (18) 

This definition is valid in the absence of magnetic inter- 
actions. Since the polarization operator commutes with 
the electrical displacement field the evolution of P(r) is 
driven only by the electronic part of Hamiltonian (jl7l) . 
The expression of the total current operator in the PZW 
gauge is: 

= :^[*nr)Vr#(r) - Vr*t(r)*(r)] (19) 
2m* 

Note that the paramagnetic term, proportional to the 
vector potential A, is absent, since in the PZW gauge 
the momentum of the particles is expressed as a func- 
tion of their velocity only [l^] . Here we have introduced 
the field operator *I'(r), constructed from the one-particle 
wavefunctions P^ : 

#(z,r||) - ^CAk</'A(^)exp(zkr||)/y^ (20) 

Ak 

Since we are interested in the intersubband transitions 
only, we shall consider solely the z-component of the cur- 
rent. The latter is readily expressed as: 

^■^^')^2S^ ^ CA,(z)e''''-|i[SA,q-<_J (21) 

A>^,q 

with the following definitions: 

<q = E^WMk (22) 

k 

Cam (2) = (t>\(z)d:,(t)^,{z) - (j)f,{z)d^(t>x(z) (23) 

For simplicity, we have chosen the envelope wavefunc- 
tions of the bound states 0a (2) to be real. Note that the 
intrasubband contribution all vanish since S,\\{z) = ac- 
cording to the above definition. 

In order to establish the polarization operator in the 
long wavelength limit, we first compute the commutators 
of the B-operators (|22p with the Hamiltonian (IT7|) : 

[^Lq' -^e] = - XI ^('^Ak+q - '^Mk)4k+qCMk (24) 
k 

In the long wavelength limit the excitation wavevector 
q is small compared to the typical electron wavevectors 
k. Since we assumed parabolic bands, we can write: 

^^Ak+q - ^A'k ~ t^A - = i^Xn (25) 
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Then iMl) becomes: 



A/jq 



Hf} — -~huj\,,B 



(26) 



This commutation relation indicates that the polariza- 
tion density operator satisfying ((T5)) is: 



Mr) = 



he 



2Sr 



E 

A>p,q 



Xfi- 



SA,,q] (27) 



This is the expression to be used in the interaction 
Hamiltonian Hint ©. It is clear that this expression 
satisfies the requirement stated in the beginning of this 
section. Indeed, it is a local function of space through 
the the microscopic current density ^^^(z). The size of 
the confining potential enters only implicitly, through the 
wavefunctions that construct £,\fj,{z) and we do not need 
to define any arbitrary averaging volume. 

To express the interaction Hamiltonian Hint as a func- 
tion of the electronic polarization, it is convenient to split 
it into two parts: 



^^int — Hii + H] 



12 



(28) 



+ 00 



ixt,iz)^X' f,' {z)dz 



(32) 



Note that the expressions derived above are exact and 
general, except for the long-wavelength approximation 
that is satisfied for the majority of experiments with in- 
tersubband devices. For instance, they will apply for a 
heterostructure featuring population inversion. In the 
following, we are interested in the specific case of non- 
inverted (thermalized) subbands. 



III. BOSONIZED PLASMA HAMILTONIAN 

A. Bright and Dark states 

In order to further study the interaction between 
the light and the thermalized subbands, we replace the 
fermionic Hamiltonian H,, by an effective bosonic Hamil- 
tonian which contains only the polarization degrees of 
freedom. In this way, we obtain a fully diagonalizable 
Hopfield-like model. To carry out this approach we need 
to replace the B-operators by effective bosonic op- 
erators. For this purpose we compute the commutator: 



with Hii the linear part and H12 the part quadratic 
in the polarization. The expression of Hn in the long 
wavelength limit, for the case of the TMq mode, is readily 
obtained with the help of dTTI) : 



(cj^kC/xk - C^k+qCAk-^q) ^ Nf, - Nx 



(33) 



Hn=i 



A>/i,q 



hu), 



cq 



-ezxfj. 



2eeoSL^ 

K-«-q)(SL-q + 5Apq) 



(29) 



Here we have introduced the dipole matrix element 
zxfj. — (0a|-z|?!'/j) of the transition /x — A, and made use 
of the following identity: 



+00 o * 

^Xp.(z)dz = zxf. 



(30) 



This identity guarantees the equivalence between the 
current and dipole matrix elements, which are more com- 
monly used in studies of intersubband transitions [l9l . [25j . 

The quadratic term of the interaction Hamiltonian be- 
comes: 



H 



12 



SeeoSm*^ ^ , ^^a^^^av' 

X>fi,X'>fi',q ^ 



E 



(<q + 5A,-q)(Bl,,,,_q + SAV'q) (31) 

where Ix^.x' ii' denotes the current-current overlap in- 
tegral: 



We recognize the difference between the number op- 
erators and N\ of the respective subbands. Let us 
define normalized operators through the relation: 



Xfj.q 



with: 



AA^Am = (Np) - (Nx) 



(34) 



(35) 



Here the mean values of the operators are computed 
in the, say, thermal state arising from the Fermi-Dirac 
statistics. In the limit of weakly excited system, the op- 
eratorsfet defined above obey the bosonic commutation 
rules 



(36) 



Taking into account ([Ml) and we replace the 
fermionic Hamiltonian He (|17l) by an effective bosonic 
Hamiltonian which yields exactly the same time evolution 
of the weakly excited system: 



E 

A>/j,q 



^^Ap&l^q&Apq 



(37) 
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FIG. 3: a) One particle excitation between two subbands. 
The fundamental subband is filled with electrons at T = 
OK. b) Bright excitation. 

The bosonic operators introduced here describe the 
only intersubband excitations that couple with light [1], 
and they are therefore called "bright" states. To illus- 
trate the bright states, we consider the system to be in 
the electric quantum limit at T = K, where only the 
first subband (A = 1) is occupied by N electrons. The 
fundamental state is then given by \F) = n|k|<fcpc|j^|0), 
where kp is the Fermi wavevector. The lowest energy 
single particle excitations can be spanned on the basis of 
states: 



|i?k> = 4kCik|F) (38) 

These basis states are eigenstates of the electronic 
Hamiltonian with an eigenenergies fvjj2i + i^fond, 
where Eiond is the total energy of the ground state. In 
the following, the energy scale is reset so that i^fond = 0. 
The excitations \E\^ are illustrated in Figure [3l^a). For 
simplicity, in our example we consider only vertical tran- 
sitions q = 0. The dipole moment operator d between 
the subbands 1 and 2 is: 



Here a,h = I..N label the new basis states. The coeffi- 
cients f3a^k are arbitrary within the normalization condi- 
tion (HIT) , but let us define the a — 1 state so that all /3i_k 
are equal. Because of (|4T|l we then have /3i_k = ^/VN 
and therefore \Bi) is exactly the state created by the 
bosonic operator introduced above: \Bi) = bl^^Q\F). 
This state, which is a coherent superposition of all pos- 
sible single particle states with equal amplitudes 1/Vn 
is illustrated in Figure |3] (b) . Next, because of (|42|) and 
(|^T|) we have: 

{F\d\Bi) = zi2Va^, {F\d\Ba^i) - (43) 

These formulas express the fact that the whole oscil- 
lator strength of the system is hold by the single state 
^q=ol-^)' whereas all other A^ — 1 states, that are orthog- 
onal to it, have a zero oscillator strength and do not 
couple with light. 

At a first glance it might seem that the emergence of 
the superradiant state 6q^Q|_F) is due to a mere change 
of the basis, and that this state does not have any par- 
ticular physical meaning. Indeed, this state is perfectly 
degenerate with the dark states as respect to the elec- 
tronic Hamiltonian Hg. However, the inclusion of the 
interaction Hamiltonian Q lifts the degeneracy between 
the bright superradiant state and the dark states \Ba^i), 
since both the quadratic H12 and the linear Hn parts 
renormalize the energy of the bright state. The quadratic 
part H12 leads to a blue shift, as explained in the next 
section, and the linear part to the emergence of two po- 
lariton states. Neither of Hji or H12 acts on the dark 
states, which therefore remain at the energy of the bare 
intersubband transition huj2i. These heavily degenerated 
dark states hinder the efficiency of the electronic injection 
into the polariton states p7| . 



= 2;i2^(c2k'^ik + cIkC2k) (39) 

k 

It is easy to show that each excitation |i?k) holds a 
dipole Z12 — {F\d\Ek). We can, however, change the 
basis by using any appropriate linear combinations of 
|^k)'s: 

\Ba)=Y,PaM\E^) (40) 
k 



E'^-A = ^a,h (41) 
k 



B. Plasma Hamiltonian 

The procedure of bosonization described in the pre- 
vious section assigns bosonic operators ^A^iq and b\^^ to 
each intersubband transition /i — > A. Each transition now 
enters the bosonized Hamiltonian with an energy huj\^. 
Therefore, in order to avoid cumbersome notations, from 
now on we label each transition /i — >■ A with a single 
Greek index a, i.e. a = [A,/i]. This means that now we 
count the number of excitations in the system, instead of 
the number of subband states. 

We now seek to express the interaction Hamiltonian as 
a function of the bosonic operators baq and 6jjq. To this 
end, we introduce the plasma frequencies tupa through 
the formula: 



{F\d\Ba) ^ Z12J2 (3a,l 
k 



(42) 



eeom*SL^ff 



(44) 
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Here L"„is the effective length introduced by Vinter 
and Tsui Q- This length is a function on the current- 
current correlation function introduced by equation (|32| , 
and describes the spatial extension of the microscopic 
current density of the intersubband transition a: 



T a 



2m*ujn 1 



(45) 



We will also make use of the transition oscillator 
strength: 



2m*Lij, 



a 2 



(46) 



The linear part of the interaction Hamiltonian then 
becomes: 



(aq-a-q)(^Lq + ^«q) 



(47) 



Here is the overlap factor between the cavity mode 
and the current distribution of the transition a, defined 
as: 



(48) 



This definition, introduced here for simplicity for the 
special case of a TMq mode can be generalized for a mode 
with an arbitrary shape, as shown in the end of this sec- 
tion. 

We now turn to the quadratic part H12 of the interac- 
tion Hamiltonian. It is expressed as a sum over pairs of 
transitions X > fj,, X' > fj,' = a, (3. Let us consider first 
the terms that correspond to the same transition, a — /3. 
Combining these terms with the electronic Hamiltonian 
([57)1 we have: 



a.q 



AuJa 



-(&iq + ^'o-q)(&Lq + &aq)] 



(49) 



This quadratic Hamiltonian can be diagonalized with 
the Bogoliubov transformation by introducing new 
bosonic operators which satisfy: 



[Paq, + -^72(0: = P)] = fu^aPac 



(50) 



where denotes the new eigenvalues. This diagonal- 
ization procedure yields the following results: 



(51) 



^ „ / ~ "aq + 1 ^ "a-q 



(52) 



In equation (jSip the new eigenvalue is exactly the 
frequency of the collective mode of the bi-dimensional 
electron gas known as the "intersubband plasmon" [2^ . 
Using (|44p and (|52|) to express the remaining a ^ /3 terms 
of Hii we arrive at the full Hamiltonian, which is now 
expressed in terms of the collective plasmonic operators: 

if = ^ ^aPiqPaq + ^ /lWcq(ajjaq -|- 1/2) 



a,q 



+ ^ nS„^(pJ,q+Pa-q)(p^_q+P/Jq) (53) 

Here we have introduced the light-matter coupling con- 
stant in the dipole gauge: 



UJPa 



L f o fw 
J a J a 



(54) 



The quantity 'E.ap is the plasmon-plasmon coupling 
constant: 



(55) 



where the coefficient Cap is the plasmon-plasmon cor- 
relation coefficient, defined as: 



\/ Ia,aIp,P \/ la, alp, P 



(56) 



The plasma Hamiltonian described in equation (|53p 
is the central result of this paper. It provides a fully 
quantum description of the coupling between the light 
and the coherent collective intersubband modes of a bi- 
dimensional electron system. The Hamiltonian (|53p con- 
tains both the interaction with the electromagnetic field 
and the coupling between plasmons from different sub- 
bands. The inter-plasmon coupling is contained in the co- 
efficients Cap- With the definition ([5^ it appears simply 
as the normalized spatial overlap between the intersub- 
band currents associated to the transitions a and /3. This 
overlap vanishes when the subbands belong to spatially 
different quantum wells. On the contrary. Cap takes val- 
ues close to unity when the subbands originate from the 
same quantum well. The coefficients Cap provide thus a 
convenient description for a number of cases, from spa- 
tially decoupled quantum wells (tight-binding approxi- 
mation), to strongly coupled heterostructures such as a 
superlattice. When several intersubband plasmons are 
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present in the system, the Bogoliubov procedure can be 
further applied to the inter-plasmon coupling terms of 
equation ()53|) in order to obtain the new normal modes 
and their coupling with the light field. These results will 
be presented in a separate paper. 

We conclude this section by providing a very general 
expression of the light-matter coupling constant for an 
arbitrary shaped guided mode fq{z)- To establish this 
expression we use the general form the z-component of 
the displacement field (O, instead of the special case of 
the TMq mode ([TT|) . We then express the linear coupling 
term Hn by following the same procedure as described 
in section IIIDI and the beginning of this section. This 
leads to the result: 



Microcavity 



/uJ^COa 6qCa^q 



(57) 




' ^12,q 



FIG. 4: Scheme of all possible interactions between a mi- 
crocavity and a three level quantum well. The thick arrow 
corresponds to the electromagnetic mode and the thin arrows 
to the intersubband plasmon modes. 



Here we have introduced the normalized current-light 
overlap coefficient: 



(|46l) and ((30|) . This leads to the expression ((54)) . where 
instead of the definition we use a generalized filling 
factor of the form: 



(58) 



and the angle 6'q is the propagation angle between the 
guided mode and the in-plane direction: 



cos Qa 



(59) 



In the previous expression e = e(z « Zq) is the back- 
ground dielectric constant of the media surrounding the 
current density (,a{z). The expressions ([55]) and ([57)1 
reveal the striking resemblance between the plasmon- 
plasmon and the plasmon-light coupling constants, S^.^ 
and Sla.q- Indeed, the two coupling constants are propor- 
tional to a normalized overlap factor, respectively Ca,ji 
and Ca^q (equations ([55t and ([551) V In the same way as 
C'a,/3 represents the overlap integral between the two cur- 
rents arising from transition a and /3, Cct^q is the over- 
lap between the current distribution a and the electri- 
cal field profile of the optical mode /q(z). Moreover, 
each plasmon enters the interaction with a weight factor 
wpa/2\/w^. The weight of the light mode is .^Wcq cos 6'q, 
the cosine term expressing the selection rule for intersub- 
band transitions |19| . Therefore, the interaction between 
the different plasmons, and the interaction with the light 
mode have the same form in the dipole gauge. Indeed, the 
Hamiltonian ([55]) describes a set of coupled oscillators, 
one of which is the electromagnetic resonator, the other 
being the collective plasmon modes. This is schematized 
in in Figure U) 

In most experimental situations the function /q(z) 
varies slowly at the scale of the intersubband current 
density ^0(2:), and therefore can be taken out from the 
integral in equation ([55)1 . In this case we recover the 
usual expression for the oscillator strength /° [2^ using 



in 



■ cos Ba 



(60) 



(The definition 



for the TMq mode is recovered 



by setting 



and Lfj 



The dipole gauge 



provides therefore a very compact description of the in- 
teraction between the light and intersubband excitations. 
In particular, the simple structure of the coupling con- 
stants allows to disclose role played by the collective ef- 
fect, described by the plasma frequencies ix)pa and the 
microcavity geometry, which sets the overlap coefficients. 
Finally, in this formulation, the weak coupling regime is 
naturally recovered for the case of vanishing plasmon- 
light overlap (Ca.q or 0), as expected. In this 
situation, which is common for absorption experiments, 
one measures solely the collective effects in the intersub- 
band system contained in the matter part of the quantum 
Hamiltonian ^ 7, 29J. 



C. Case of OD microcavities 

So far we used the expansion of the electromagnetic 
into the basis of guided modes labelled by the wavevec- 
tor q. Our choice was motivated by the microscopic def- 
inition of the polarization density in section III D[ which 
is naturally expanded into elecronic plane waves propa- 
gating along the hetcrostucture slab (see, for instance, 

m)- 

The approach developed here allows the rigorous quan- 
tum description of OD microcavities that confine the 
electromagnetic field into all tree dimensions of space. 
Such systems have been recently employed for the study 
of light-matter coupling with intersubband transitions 
@, [30} . In this case the possible polarization excitations 
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will be determined by the quantizing conditions for q 
imposed by the microcavity boundaries. 

We start by expanding the electrical displacement into 
laterally localized TMo-like modes: 



D(r) = D(z,r|| 



25Lc 



-u,„(r||)(aj'„ - a™) (61) 



The index m labels the discrete cavity modes and the 
set of functions Um{v\\) describe the lateral shape of the 
modes. They are normalized such as: 



Um(r||)wm'(r||)c^^i"|| = 'S'<5r, 



(62) 



We must now provide the expansion the polarization 
density (P7| into the new basis of localized electromag- 
netic modes. To this end, it is very convenient to use 
the quantum mechanical notations |q) and \um) so that 
(r|||q) = exp(iqr||)/-\/^ and (r|||um) =Wm(r||). Then the 
basis transformation is expressed as: 



(63) 



Here (q|um) is simply the q*^ Fourier component of 
the function Um(r||)[3l|- 



{q\um)^\ jj e^^^\^uUv\\)d\\ (64) 

We have, accordingly, the transformation law for the 
B-operators (f22|) : 



am 



(65) 



a) 



'-cav 



Electric displacement D 



b) 




QW 



Intersubband charge 
density (1 2) 



FIG. 5; a) Electrical displacement for the (iV = 1,M = 0) 
mode of a patch microcavity with a patch size s. b) In- 
tersubband charge density for the fundamental transition 
of a quantum well with thickness Lqw coupled with the 
{N = 1, M = 0) mode of the microcavity. 



The above Hamiltonian (1571) is formally similar to the 
plasma Hamiltonian (|53p . The difference arises from the 
modified spectrum of electromagnetic modes, which is 
no longer continuum but discrete. The discrete spec- 
trum allows to study the simplest possible light-matter 
interacting system, where only the lowest intersubband 
transition [a — 2,1) interacts resonantly with a single 
microcavity mode that is sufficiently far from the oth- 
ers, say the fundamental mode m — Q. Then the plasma 
Hamiltonian for this "polariton dot" system is reduced 
to [S]: 

H = /iW2iP2ii'2i + hujc{a)a + 1/2) 



+i 



- J^f2if2i{a^ -a){pl,+p2i) (68) 

V ^21 



An example of a microcavity is provided by the square 
patch resonators depicted in Refs. [l,!!^]- In this case the 
resonances are labelled by two lateral indexes m = N,M 
and the corresponding lateral functions are: 



Since the same transformation law applies to the 
bosonic operators &jjq, we can readily express the expan- 
sion of the polarization density in the new basis: 



he ^ ^^" ^a(z)M„,(r||)[b^,^ 



2Sm^ ' — ' uj, 

OL.ra 



1(66) 



The interaction Hamiltonian is also readily expressed. 
Details are provided in Appendix [B] and the final ex- 
pression for the plasma Hamiltonian of OD microcavities 



^ = ^ flOJapi,„Pam + ^ huJcm{at,am + 1/2) 
a,m m 

a.m * 



■"Ar,M(l"||) = 

^2 - SonV^ - 5oM cos (^) cos (^^) (69) 

Here s is the size of the square. The fundamental mode 
is actually twofold degenerate (iV = 1,M = 0),(A^ = 
0, M = 1), however only one of the two modes can be 
selected in experiments with polarized light. The mode 
(iV = 1,M = 0) is depicted in Figure [Sla) . 

Generally, in the case of a resonant light excitation, 
only the spatial modes that correspond to the specific mi- 
crocavity resonance will be excited. That means that the 
charge distribution will vibrate into well defined spatial 
modes, imposed by the microcavity. The charge density 
p(r) can be obtained from the divergence of the polariza- 
tion matrix element p(r) = dPz{v)/dz: 



.(r) 



A/1 



(j)x{z)(j)t,{z)uyyi{r\\) (70) 
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In the case 
of thickness i. 



of 

QW 



a square infinite quantum weU 
the wavefunctions are 4>\(z) — 
sin(A7r2:/LQw)- The charge distribution which 
corresponds to the first intersubband excitation 1 — )• 2 
coupled with the fundamental cavity mode N — 1,M ~ Q 
has been plotted in Figure[51Jb). This example shows how 
the microcavity allows to control the spatial properties 
of the bright state. 



D. Back transformation in the Coulomb gauge 

The light-matter coupling between microcavities and 
intersubband transitions have been studied, so far, ex- 
clusively in the Coulomb gauge Q . Usually, the case of a 
single intersubband transition coupled with a continuum 
of guided modes Wcq has been considered. These theo- 
retical studies pointed out the possibility to obtain the 
"ultra-strong" coupling regime, where the light- matter 
coupling constant becomes comparable to the frequency 
of the intersubband transition 0^21 [3]. In this case the 
full quantum Hamiltonian, including the anti-resonant 
terms and the quadratic vector potential term must 
be taken into account in order to describe correctly the 
system Q. 

In this section we establish a link between the previ- 
ous studies in the Coulomb gauge and the description in 
the dipole gauge developed here. We shall consider the 
case of a single intersubband transition (A = 2,/i = 1), 
coupled with the TMq mode. We then perform a unitary 
transformation to the Hamiltonian in order to obtain its 
expression in the Coulomb gage. The transformed Hamil- 
tonian will have identical eigenvalues as the original one, 
but will be expressed in terms of the vector and scalar 
potentials A and V. 

In this section, it is convenient to keep the individ- 
ual subband indexes 2, 1. We return to the form of the 
Hamiltonian with a single intersubband transition before 
the Bogoliubov transformation leading to the depolariza- 
tion shift: 



H^J2 ^2i6^foq + Yj ^cq(a^aq + 1/2) 



+ ^ E ^ - + 6q) 

fe, ,2 

+r^(^l + ^-q)(^-q + M (71) 

(The subscripts "21" have been dropped for operators). 
The vector potential for the TMq mode is: 



T = exp ( - J A(r) • P{r)d^r^ (73) 
It writes in the case of a single intersubband transition: 



T = exp ( - i ^ Xq(aq + a-q)(&-q + bq) 



^1 / ^p/21/21 
21/ a;2iWcq 



(74) 
(75) 



Recalling the parity of the photon dispersion Wcq = 
the following transformation laws of the bosonic 



^c— q 



operators are obtained: 



r+6qT = &q - iXq(a-q + Aq) 
r+flqT = aq - «Xq(6lq + foq) 



(76) 
(77) 



With these relations and their hermitian conjugates, 
the transformed Hamiltonian is obtained to be: 



T+HT = ^ huj2iblbq + fi^cq(a]jaq + 1/2) 



-i Y ftf^q(aq + a-q)(&q " ^-q) 



HQ 



E — ^("q + «-q)(«-q + oq) 



(1 - /2l/2'"l) E r7^(^q + ^-<l)(^-q + ^) (^8) 



Here Qq is the light-matter coupling constant in the 
minimal coupling gauge: 



q = '^2lXq 



UJp UJ21 



'J21J21 



(79) 



^cq 



The first four terms of (|78t , together with ([79| provide 
the Hopfield-like Hamiltonian used so for the theoret- 
ical study of intersubband polaritons d, H^. Namely, 
the presence of the A^ term and the anti-resonant terms 
leads to the " ultra-strong" coupling regime Q . The later 
has been defined as the situation where the coupling term 
57q, taken at resonance W21 = Wcq becomes comparable 
with the energy of the intersubband transition uj2i'- 



A(r) = ^ , 



2eenSL 



O'J-'Jcav'-^cq 



^''^■•Haq + aLq) (72) 



We use (TTp to express the inverse PZW unitary trans- 
formation [l|: 



nqiu;21 = UJ,q) = ^Vf2j2l ~ ^21 (80) 



Since from (ffS]) and ([79|) the A^ term is proportional 
to square of the plasma frequency Wp, we see that the 
ultra-strong coupling regime, as defined in Ref. Q is 
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obtained in systems featuring high photonic confinement 
factor f2i, and high plasma frequency ujp, i.e. high elec- 
tronic densities. However, for high electronic densities 
the Coulomb interaction bring important dynamical cor- 
rections ^8]. The latter are already present in the dipole 
gauge, which naturally includes the collective excitations 
of the electron g described in section IIII Bl In the 

Hamiltonian ([75)) expressed in the Coulomb gauge, these 
correction actually arise from the last term. In Appendix 
[Q we show that this term can indeed be cast in a form 
of a long- wavelength limit of the Coulomb potential. 

The Coulomb correction in the Hamiltonian (|75)) con- 
tains an interesting new element, which is the geometrical 
prefactor (1 — /2i/2i)- This prefactor is equal to one in 
the case of a very large cavity {f2i — >■ 0). This case of 
vanishing photon confinement corresponds, for instance, 
to the multi-pass waveguides employed for the absorp- 
tion measurements where the depolarization shift is ob- 
served. However, the confinement factor becomes 
an important correction in the case of micro-cavities with 
filling factors f^i close to unity, such as the double metal 
microcavities [6, 11, 30]. We interpret the factor — /21/21 
as an image contribution to the Coulomb interaction due 
to the boundary conditions of the electric field at the mi- 
crocavity walls. In other words, it can be seen as a local 
field correction due the partial screening of the microcav- 
ity field by the oscillating intersubband charges. Indeed, 
both the displacement field ^^(r) and the polarization 
field Pz{j^) that we used to construct the PZW Hamilto- 
nian satisfy the boundary conditions on the cavity walls, 
which are transported to the Coulomb correction through 
the unitary transformation ([74l) . 

Note that the the unitary transformation ([74l) holds 
only for a truncated Hilbert space, where we retained 
only the first intersubband transition and the fundamen- 
tal waveguide mode. The full gauge equivalence is es- 
tablished correctly only if the complete set of quantum 
transitions and electromagnetic modes of the system are 
accounted for 34] . This issue will be discussed elsewhere. 

We can combine the last term in the Hamiltonian ([78]) 
with the first one, and perform a Bogoliubov diagonal- 
ization just like in section flH Bl As a result, we obtain 
a renormalized intersubband frequency with an effective 
depolarization shift, which takes into account the local 
field corrections, contained in the factor /21/21: 

^21 = yJ^li+^Ul- f^J^i) (81) 

The local field factor 1 — yields an effective 

plasma frequency wp-^/l — /2i/2i- We can then rewrite 
the full Hamiltonian in the Coulomb gauge ([75)1 using 
equation (|8ip and the corresponding polarization opera- 
tors. The result is exactly the Hamiltonian of Ref. 
used for the study of the ultra-strong coupling regime, 
where the bare intersubband frequency has been replaced 
by the effective frequency uj21- 



T+HT = mib\b^ + fi^cq(al,aq + 1/2) 
q q 

q 

+ E^K + «-q)(«-q + «q) (82) 

Note that now 6q describes the bosonic operator after 
the Bogoliubov transformation, and the frequency 0)21 
should be used in the definition of Hq. This result vali- 
dates the studies performed in the Coulomb gauge [ssl ] , in 
the limit of microcavities with large filling factors, where 
W21 ~ UJ21- However, contrary to Refs. [1] and [ssj . 
the effective Hamiltonian (|82l) includes also the correct 
limit of the weak coupling regime, obtained for vanish- 
ing overlap with the cavity mode — >■ 0. In this case 
from expression (|8ip we recover the renormalized tran- 
sition frequency CD21, whereas the Hamiltonians in Refs. 
Q and would predict only the bare intersubband 
spacing ti;2i- 



IV. PROPERTIES OF THE POLARITON 
STATES 

A. Polariton dispersion 

We now analyse coupled light-matter polariton states 
arising from the dipolar Hamiltonian. To simplify, we 
ignore the coupling between plasmons on different sub- 
bands. This is equivalent to consider a single intersub- 
band transition uja, not necessarily the fundamental one, 
in interaction with the TMq mode. All the results that 
will be stated remain also valid for OD resonators. The 
corresponding plasma Hamiltonian is: 

H = fl^aPqPq + ^ ^^^cqlaqAq + 1/2) 

q q 

+*E^\/S^(at-a-q)(pLq+Pq) (83) 
^ z y 

This Hamiltonian is very similar to a Dicke model 
[36j . However, the coupling coefficient is proportional 
to Ldp/\/u^ and has a non- linear dependence on ujp be- 
cause of the formula of the depolarization shift uja — 
+ ■ We show further that this non-linearity leads 
to the no-go theorem for quantum well systems, and the 
Hamiltonian (|83p does not allow a quantum phase tran- 
sition ^]. The Hamiltonian (j83| can be diagonalized 
exactly by introducing the polariton operator: 

Hq = Xqflq + J/qflLq + ^qPq -|- tqpLq (84) 
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The Hopfield coefBcients introduced here satisfy the 
normaUzation condition: 



\x^\'-\yc\' + \z^\'-K\' = i 



(85) 



The Hopficld-Bogohubov determinant corresponding 
to the equation [H, Ilq] — fiwqllq is then: 



cq ~ '^q 







if2q 





—'^cq — l^q 


ifiq 


ifiq 






" Wq 







— if2q 








(86) 



with ilq the hght-matter couphng constant from ([M)) 
(the subscript a has been dropped). Zeroing the deter- 
minant ([86)1 provides the eigenvalue equation: 



cq/ 



2 , 



cq 



(87) 



This biquadratic equation can be solved analytically, 
the two real solutions, u;q,+ and a;q^_, being the frequen- 
cies of the two coupled states. The Hopfield coefficients 
are also readily expressed in closed form. For instance. 



we can define a "photonic" 
an "electronic" part — \ 



part hp 



and 



|fq| linked by the re- 



lation hp + I 
expressions: 



hp,+ — 



1. For the photonic part we obtain the 



Note that we have necessarily hp^. 



therefore he + + h, 



hp- 



(88) 



1 (and 



The polariton frequencies ojq^±, as well as the elec- 
tronic Hopfield coefficients he.± have been plotted as a 
fimction of the cavity frequency Wcq on Figure [6]^ a) (c). 
For this illustration we have used the numerical values 
Ljp = 0.83w„, /° = 1 and = 0.8. As seen from Figure 
EUa) the polarion dispersion features a gap. The upper 
edge of the gap, obtained at Wcq = is the frequency of 
the intersubband plasmon uja- The lower edge is easily 
estimated from equation ([57)1 to be: 



fo fw\ 



(89) 
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FIG. 6: (a) Polariton dispersion, normalized at the bare in- 
tersubband transition lOa (b) Splitting of the two polariton 
states, (c) Hopfield coefficients. In these figures, the sub- 
script q has been dropped, for clarity. The numerical values 
used for solving equation H87p are top = O.SSoJq., /a = 1 and 
/a = 0.8. 



and photonic degrees of freedom. This mixing is quanti- 
fied by the Hopfield coefficients (|88)) . Namely, the cou- 
pled system features maximum mixing when the photonic 
part of the Hopfield coefficients equals the electronic part: 
^e,± = hp_± = 1/2. From equation ([55)1 we readily ob- 
tain that this is satisfied when the cavity is tuned into 
resonance with the intersubband plasmon: 



hp,± ^ t^cq 



UJn 



(90) 



This has been illustrated in Figure Efc). In the well 
known Jaynes-Cummings model |37| the maximum mix- 
ing also corresponds to the point of minimal splitting 
a;q_+ — Wq.- between the two polariton states. However, 
this is not true in the general case. As shown in Ap- 
pendix [Dl the minimum splitting occurs when the cavity 
is resonant with the lower gap edge frequency (equation 
dHU)): 



In the next section we show that /°/™ < 1, therefore 
the lower gap edge always appears at a frequency higher 
than the bare intersubband frequency uja- The impossi- 
bility for the light to propagate at the gap energies can be 
explained by the destructive interference between the mi- 
crocavity electromagnetic field and the local field created 
by the collective electronic oscillations. 

The most important aspect of the strong-light mat- 
ter coupling regime is the mixing between the electronic 



d(a;, 



q,+ 



duj. 



= LUcq = ijJa 



(91) 



cq 



This is also illustrated in Figure Hfb). Moreover, in 
Appendix |D] we show that the minimal splitting can be 
computed exactly, and the results is: 



min(a;q.. 



-) - V^^P = 2r!fi (92) 




FIG. 7: Polariton frequencies as a function of the plasma 
frequency for a cavity resonant with the bare intersubband 
transition ojcq = uja- The grey area indicates the gap. For 
this figure we have used fa = 1 and /™ = 0.8. The linear 
asymptotes correspond to the Jaynes-Cummings model. 

We recognize the quantity 2^2 to be the Rabi splitting, 
as known from the usual definition f38|: 

Note that this is an exact result, that is valid for all 
orders of light- matter interaction. More general result for 
an arbitrary cavity can be obtained with the expressions 
([57)1 or ([50]) . The quantity 2ilfj can therefore be used 
as an experimental measure of the interaction strength, 
even in the ultra-strong coupling regime. 

We see that both edges of the gap, <Dq and oja play an 
important role in the theory. The frequency uja defines 
the point of maximum quantum mixing, whereas Wq. de- 
fines the point of minimum splitting 2f2ij. However, we 
can observe from the plot on Figure [6ljb) that the split- 
ting around the point Wcq = is not very different from 
251^. This feature is a characteristic of the ultra-strong 
coupling regime, where we can no longer define a strict 
resonance condition for the optimal coupling point. In- 
deed, the system will feature almost identical coupling 
energy 2Mlji if the cavity resonates with any frequency 
in the band between uia and Wq . The three characteristic 
frequencies w^, uja and 2VLji are liked through the simple 
relation: 

Z:l=u:l+ml (94) 

As pointed out in Ref. for a very large interaction 
strength, the effects of the quadratic and anti-resonant 
terms in the interaction Hamiltonian become important. 
These effects are manifested with the non-linear behav- 
ior of the polariton frequencies Wq_-|- as a function of the 
Rabi splitting 20, n [4]. Since the light-matter interaction 
strength scales as the plasma frequency ujp (equations 
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FIG. 8: Illustration of the polariton frequencies for cavities 
with increasing filling factor /™. For this illustration, the 
cavity is slightly red-shifted as respect to the bare cavity fre- 
quency LJa- 

((92)) and l|93|) ') we have studied these effects as a function 
of Up, as shown in Figure [71 Experimentally, ujp can be 
varied either through the temperature of the system, or 
by applying a gate voltage [39|. In both cases one con- 
trols the subband population difference ANa- For the 
plot of Figure [7] we have chosen a cavity that is resonant 
with the bare intersubband transition when ujp = 0, i.e. 
i^cq = Wq. Moreover, in this example, we neglect the 
Hartree correction of the bare intersubband frequency as 
the number of charges in the system progressively in- 
creases. This assumption is true for a sufficiently thin 
square quantum well. Because both frequencies Wq. and 
uJa increase with ujp due to the depolarization effect, 
the system gets blue-shifted from the cavity mode, how- 
ever the blue shift is much smaller for uja- The Jaynes- 
Cummings model corresponds to the linear asymptotes 
of the polariton branches at low utp. Figure [7| clearly 
shows that this model ceases to be valid for high elec- 
tronic densities. Moreover, the magnitude of the polari- 
ton gap increases as the polariton branches depart from 
the linear asymptotes. Therefore the measurement of the 
gap can be considered as a direct spectroscopic evidence 
of the ultra-strong coupling regime [ll| . 

Along with the charge density, the second ingredient 
of the light -matter interaction is the overlap factor /™ 
between the photonic mode and the intersubband cur- 
rent. As evident from the expression of the Hamiltonian 
in the Coulomb gauge ([75)) . this factor controls not only 
the intensity of the light-mater coupling (linear term in 
equation ([75)) ). but also the relative weigh between the 
longitudinal Coulomb corrections (last term in ([75]) ) and 
the transverse corrections contained in the term. 

The influence of the overlap factor on the polari- 
ton dispersion is illustrated in Figure [5| In this plot, the 
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FIG. 9: Sketch of the intersubband current £,c,{z) for a het- 
erostructure in a cavity witli thickness Lcav ■ We suppose that 
the current decays sufficiently fast, so it can be considered as 
null close to the cavity boundaries. We can then define a re- 
gion of space with a thickness Lqw < icav containing the 
entire current density. 



cavity frequency ojcq is slightly blue-shifted as respect to 
the bare intersubband transition Wq. When /™ = 0, we 
recover two uncoupled oscillators at frequencies ujcq and 
uja as expected from equation (|571) . For small values of 
there is a small splitting that appears when the cav- 
ity is tuned with the intersubband plasmon, w^q = uJa- 
In this case the term of the Hamiltonian ([7§|) . that 
scales as )^ is negligible. The polariton gap is also 
negligible, since lDq, ~ Wq, and the system is described by 
Jaynes-Cummings model with a resonant frequency ujq, 
renormalized by the Coulomb interactions. 

On the contrary, for filling factors close to one, ~ 1 
the Coulomb correction in the Hamiltonian ([75]) is neg- 
ligible, and the weight of the quadratic term shifts to 
A^. If the electronic density is sufficiently high, the sys- 
tem enters the ultra-strong coupling regime as defined in 
Ref. 3- 

The resonance condition Wcq = looses its 
strict meaning, as the cavity can be resonant with any 
frequency in the polariton gap. 



B. No-go theorem 

The "No-go" theorem states the impossibility of the 
lower polariton state to acquire zero energy when the 
light- matter interaction is increased [T2j . This property 
is related from the quadratic term in the interaction light- 
matter Hamiltonian. For intersubband polaritons, the 
No-go theorem becomes particulary clear in the dipole 
gauge, where, as shown bellow, it stems from the depo- 
larization effect. 

The asymptotic value of the lower polariton frequency 
for very high plasma frequencies — )• oo is deduced 
from equation (|87p : 



(95) 



Note that this limit does not contain the bare intersub- 
band frequency uja , and therefore is independent from the 
eventual Hartree corrections to the heterostructure po- 



tential. The No-go theorem for intersubband polariton is 
then equivalent to the following strict inequality: 



J a J a ^ 



(96) 



This inequality can be proven by using the properties 
of the intersubband current matrix element (,a{z) 
which allows to express the overlap factor /™ — L^^/Lc 
through equation ([iS)) : 



rw 



(97) 



We suppose that wavefunctions of the bound states 
and their derivatives decay sufficiently fast, so that fa (2) 
is zero close to the cavity boundaries, as illustrated in 
Figure IHl Then we can define a domain with a finite size 
^QW < icav, that we can arbitrarily call "quantum well" , 
such as all the the wavefunctions and their derivatives 
are zero outside this domain (Figure [2). We chose the 
origin of the coordinates so that ^ z ^ ^qw inside the 
" quantum well" . We then re-express the current-current 
integral in the non-local form: 



+00 



+ 00 



(98) 



We can choose an arbitrary orthogonal basis of real 
functions Xn{z) on the segment [0, Lqw] in order to span 
the delta-function: 



5{Z - 2') = X! ^^(^)Xn{z'), 



(99) 



This basis is not necessarily the basis of envelope wave- 
functions. The expansion of the integral (|98| on the basis 
Xn{z) is: 



/ U^fdz^J2{ Uz)Xn{z)dz) (100) 

An evident choice for Xn{z) is the Fourier basis on the 
segment [0, Lqw]: 



Xniz) 



Lqvj 



Lq-w 



■ sm 



( 2-!Tnz \ 
Lq-w J 



odd 



(101) 



With this particular choice, the first basis function 
Xo{z) = -y/l/Lqw is constant on the segment [0,Lqw]. 
Since all terms of equation (jlOOp are positive, we have 
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FIG. 10: Comparison of the polariton dispersion, at resonance 
(ucq — uja), with and without the depolarization effect. For 
simpUcity, the Hartree shift of LUa is neglected in this illustra- 
tion. 



the Cauchy-Schwartz inequality for the current-current 
overlap: 



-t-oo 



Uz)dz) (102) 



The right hand side of the last equation can be ex- 
pressed with the oscillator strength according to the iden- 
tities (pO)) and (|46p . Combining this expression with the 
definition (j97p . we obtain the following inequality: 



Lqw 
f° 

J a 



(103) 



This result shows that the higher the oscillator 
strength of an intersubband transition, the more the cor- 
responding current density is localized in the space. This 
provides directly the No-go theorem, since it leads to the 
inequality: 



/w ro ^ 
a J a 



L 



QW 



< 1 



(104) 



Therefore the No-go theorem is a consequence of the 
confinement of the electronic plasma in a finite volume 
of space [l^l . This confinement is also the origin of the 
depolarization shift, which arises from the plasma energy 
of the bound charges. 

The inequality (jl03p provides a lower bound for the 
plasma energy of the system. The plasma energy is quan- 
tified as the weight of the term in the plasma Hamil- 
tonian (gH]): 



fe2 



ASL 



(105) 



Using the inequality (jlOSI) and the definition of the os- 
cillator strength (|46|) . we obtain the following inequality 
for the plasma energy: 



> 



2eeo 5"^ 



QW 



(106) 



Here da = &Za is the dipole moment of the inter- 
subband transition. The left side of (|106p corresponds 
exactly to the classical self-interaction energy of AiVo, 
dipoles uniformly distributed in a volume 5'Lqw [13 ■ 
According to equation ()106|) the quantum self-energy is 
higher than the classical one, since the wavefunctions of 
the bound states illustrated in Figure [2] provide stronger 
spatial confinement of the charged particles. 

The plasma energy of the charged particles localized 
in a finite volume compensates for the energy decrease of 
the lower polariton as the light-matter coupling strength 
is increased. This is illustrated in Figure [TUl where we 
have compared the polariton branches w± as a function 
of LUp with the solution of equation (1571) when the de- 
polarization shift is discarded, i.e. when uja is replaced 
with LUa ■ The cavity is chosen so that Wcq = Wq , and the 
Hartree corrections have been neglected for clarity. It is 
clearly seen that the depolarization effect prevents the 
lower polariton energy to reach zero, which is an equiv- 
alent statement of the No-go theorem for intersubband 
transitions. 



Effective dielectric function 



The dispersion relation ([87| allows to obtain the effec- 
tive dielectric constant ecs{ui) of the polaritonic medium. 
The latter can be defined through the propagation equa- 
tion: 



eoff(wq)-^ = 



(107) 



Since for the TMq mode cj^^ = c q /e, we obtain from 
83: 



ecff(w) 



= 1 



fo fw 2 



(108) 



Note that the zero of the dielectric constant Ecff (<^) cor- 
responds to the frequency of the intersubband plasmon 
Wq, as expected. The above expression can be recast in 
the form: 



1 J^O -fv 

J a J a 



eoff(w) 



£Qw(t^) 



(109) 



off 



where £qw('^) is the dielectric constant of the het- 
erostructure alone: 
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CONCLUSION 



eQwM^.(l-^) 



(110) 



The result p09p is very similar to the effective dielec- 
tric constant obtained by Zeluzny et Nalewajko [3] ■ This 
result arises naturally from the initial assumption that 
the heterostructure is much smaller than the wavelength 
of light. However, in our case the geometrical overlap 
factor f°f^' contains a quantum correction due to the 
shape of the wavefunctions. 

The model developed so far is purely Hamiltonian, and 
the dissipation has not been included. The dissipation 
can be introduced through a coupling with a reservoir of 
harmonic oscillators. Such model has been considered, 
for instance, in Refs. [13, US- In particular, the method 
of Ref. UM leads to a dielectric function of the form: 



£qw(w) = £( 1 



Pa 



(111) 



a;2 + iTaUj — uj"^. 

Here Ta is the linewidth of the intcrsubband plasmon 
modes, that can be determined in a phenomenological 
way, i.e. from absorption measurements |19| |. 

On the basis on the results from the previous section we 
can obtain an approximation for the dielectric constant, 
by replacing the effective plasma thickness _L"g with 
its maximum value iqw/Za (see equation (|103p ). We 
can call this substitution " semiclassical approximation" , 
since it leads to a plasma self-energy provided by the clas- 
sical expression, as seen from equation ()106|) . This leads 
to the widespread expression for the heterostructure di- 
electric constant [23: 



(112) 
(113) 



£QW,ci(w) = £ 1 - 



eeom* SLqw 



In this limit, we obtain exactly the effective medium 
constant from Ref. [l^ . since in this case the product 
fafa is replaced with iqw/^^^cav, which is independent 
from the particular intcrsubband transition. From sec- 
tion |IVB] it is clear that the semiclassical approximation 
assumes that the intcrsubband polarization is constant 
inside the heterostructure slab, and zero everywhere else. 

Including more than one transition in the definition 
of the dielectric constant is not not trivial, since one is 
obliged to diagonalize the full plasma Hamiltonian ([55)) 
taking into account the coupling between the different 
intcrsubband plasmons. The semiclassical approxima- 
tion allows to write a simple analytical expression for 
the mutli-transition dielectric function: 



eQW,ci(, 



■"Pa 



This result is demonstrated in Appendix lEl 



(114) 



In summary, we have provided a theoretical description 
of the light-matter interaction for the intcrsubband tran- 
sitions in the electrical dipole gauge. We showed that, 
by introducing a microscopic expression for the intcrsub- 
band polarization field P(r), the Power- Zienau-Woolley 
Hamiltonian provides a suitable framework for the study 
of the interaction between the collective electronic exci- 
tations and light. This description is very general and 
applies to an arbitrary electronic potential, once the cor- 
responding single particle wavefunctions are known. In 
particular, it can be used in the case of the ultra-strong 
coupling regime, in the limit of very high electronic den- 
sities and high overlap factors between the quantum well 
medium and the microcavity mode. 

The physical interpretation of the dipolar interaction 
Hamiltonian is straightforward. The linear part of the 
Hamiltonian, Hn describes the coupling of the electronic 
excitations with the cavity electromagnetic field. The 
quadratic part, -ff/2, contains the depolarization effect 
(the effect of an oscillating current on itself) and the cou- 
pling between plasmons from different transitions. The 
plasma energy described by the term H12 leads to the 
No-go theorem for intcrsubband transitions. This energy 
has a close analogy the electromagnetic dipole self-energy 
of a two level system [lo| . that plays an important role 
for the estimation of the radiative corrections of atomic 
transitions [l|, . 

When the dipolar Hamiltonian is transformed back in 
the Coulomb gauge, it provides the Coulomb interaction 
terms of the system including image charge effects arising 
from the boundary conditions on the cavity walls. This 
has been discussed for a truncated Hamiltonian describ- 
ing a single intcrsubband transition interacting with a 
TMq mode of a double metal waveguide. In this case we 
find that the weights of the scalar potential V and the 
term of the Coulomb gauge version of the Hamiltonian 
are governed by the geometrical overlap factor /™. Our 
study completes the theoretical framework of the ultra- 
strong coupling, introduced in Ref. 4], where only on 
the vector potential part A^ was considered. We found 
that this study is justified for cavities with filling factors 
close to unity. A more general PZW back transforma- 
tion, which uses the whole basis of guided modes and 
the full set of intcrsubband transitions will be discussed 
elsewhere. 

We believe that the approach developed here, based 
on operator formalism to describe the collective intcrsub- 
band excitations, provides a compact and useful descrip- 
tion that could allow further studies of quantum electro- 
dynamical effects inserted in solid state systems. 

The authors acknowledge very useful discussions with 
C. Ciuti, D. Hagenmuller and S. De Liberato. We also 
acknowledge financial support from the ERC grant " AD- 
EQUATE" . 



18 



Appendix A: Vacuum field amplitude of the guided 
modes 



In this Appendix, we derive the vacuum normahzation 
constant Aq (equation ([7])) by taking into account the 
multilayered geometry of the guided mode stack. We also 
comment on the special treatment of metallic boundaries 
of the guiding structure. 

We start by expressing the volume integral ([2]) as a 
function of the TM polarized field components (4)-(6): 



Here we have also used that eoMo = and the 

symbol i± means the upper/lower side of the layer. Using 
the boundary conditions of the electromagnetic field: 



1^ 

£,■ dz 



/q|i- — /ql(i+l) 
1 rf/q 



dz 



(■'+!) + 



(A7) 
(AS) 



and the fact that 7=1/2 J^ - finally obtain: 



Mo 
2 



S ^ Aq(a]jaq + aqolj + Aqfl-q + aqflLq) X 

q 

^ j fl{z)dz (Al) 



-D?d^r 



2eae{z) 



2eoe{z)' 

S ^ ^q(aqaq + OqC^ - aqfl^q - ajjoLq) X 



1 



^cq J 2£oe(z) 



fl{z)dz (A2) 



2eae{z) II 



S ^ ylq(a]jaq + a^al^ - a^a^^ - a\a}_^ x 



1 /rf/q 
cj^ J 2eo£{z) V dz 



dz (A3) 



To derive these expressions we have used the parity of 



^cq 



and Ac^ as respect to q and e_ 



-Bq. We have 



also used the orthogonality of the two dimensional space 
harmonics exp (iqr|| ) / ^/S. 

To obtain the proper normalization of the electromag- 
netic field, we first show the following identity: 



dz 



2 ■''^ 2eoe{z) \uj, 



dz 



(A4) 



For the demonstration, we use the Helmholtz equa- 
tion Let us consider the i^^ layer, then by definition 
e{z) — £i is constant. By multiplying the equation (|9]) 
by /q and integrating across the i"^ layer we obtain: 



dz 



S{)Si dz 



Mo/q 



1 



i+ 



— f 



S()£i ^ dz 



_f ^ 



dz ) 



(A5) 
(A6) 



1 



2eaeu dz 



1 J dh 

2e()SD dz 



(A9) 



where U/D denote the upper/lower boundary of the 
multi-layer. Either one or the two boundaries go to in- 
finity; then the field and its derivative go to zero. The 
other option is to take metallic boundaries. Since in the 
mid- and far- infrared range the dielectric constant of 
the metals Em is very high, we can adopt perfect metal 
boundary conditions, where {emI oo. We can then ne- 
glect the contribution of the field from the metallic layers 
in the integral ([2]). In all cases the identity (jA4| is proven. 
Thanks to this identity, the anti-resonant terms in the 
photon Hamiltonian cancel, and the remaining term is 
simplified to: 



q 



2^o^iq^q(a]jaq + 1/2) 



(AlO) 



Following standard textbooks [10|, the prefactor 
should correspond to the energy quantum of the guided 
mode hujcq, which leads to equation ([7|). 

Although we used perfect metal boundary conditions 
for the normalization of the vacuum field amplitude Aq, 
we can still describe the effects of the finite metal per- 
mittivity £m in our model. As an example, consider the 
current experimental situation where the multilayer is 
constituted by a surface plasmon waveguide, obtained 
by deposing a metal layer on a semiconductor substrate. 
We suppose the metal-semiconductor interface to be at 
z = 0. In the semiconductor the mode profile is described 
by an exponential function fq{z) — e'^^^ for z < 0, with 
7q the decay length of the surface plasmon mode into 
the semiconductor. Let e be the dielectric constant of 
the semiconductor, then we have |4l| : 



2 

cq 

,2 



= q 



7q 



cq 



(All) 
(A12) 



Using the above prescription, we consider the field in 
the metal to vanish, but we keep the expression of the 
finite dielectric function Em in the dispersion relations 
above. We then have Lq = l/27q, and the quantum 
amplitude of the surface plasmon mode becomes: 
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1/2 



(A13) 



We see that in the far infrared domain, where — >■ 
and Em is very large and negative, both spontaneous 
emission and the strong couphng phenomena are un- 
favoured, since has vanishingly small values. On the 
contrary, the photon electric amplitude can be en- 
hanced in a microcavity resonator. 



Appendix B: Derivation of the polartion dot 
Hamiltonian 



We first write the linear part of the interaction Hamil- 
tonian ^ by using the expressions (P?]). (pT|) and the 
definition (l6il): 



hljJr 



-Za{cily^ dm) X 



2eeoSL^ 

UL,ftL 

Y,[Bi^{^\Urn) + Baq(-q|Wm)] (Bl) 

q 

Since the function Um(r||) is real, we have (— q|itm) = 
(wm|q), and therefore we can use the transformation law 
(|65p and its hermitian conjugate to obtain: 



huJr 



2eeoSL^ 



{Bam + Bo 



(B2) 



In order to transform the quadratic part (j3ip we use 
the inverse transformations: 



BU = T.^^m\ci)B^ 

m 

Baq = y^^{q\um)B, 



t 

otm 



(B3) 
(B4) 



and therefore: 



Y,iBU + Ba-^)iBl_^ + Bp^) = 
q 

{{um\q}Bl^ + {-q\Um)Bcm) X 

q,m,m' 

((w„,'|-q)5^„, + (q|u™')5/3m') (B5) 

Using once again the relation (— q|u„i) — (u,„|q) and 
the closure relation J2q |q)(q| = ^ along with the orthog- 
onality condition {um\um') = ^mm' it is straightforward 
to transform the expression (|B5I) into: 



Y^{Bl^ + Bam){Bl„,+B, 



Pm) 



(B6) 



We then use exactly the same bosonization approach 
as in section UlI Al to arrive at the Hamiltonian (|67|) . 



Appendix C: Long wavelength limit of the Coulomb 
interaction 



In order to show that the last term of equation (|78l) has 
the form of the Coulomb interaction, we first re-express 
it through the _B-operators defined in section Hi D I 



5: S'^i 



2t2 



■/l2,12 E(^21q + S21-q)(Sll_q + i?21q) (CI) 



We first use the Vinter's indenity [tI, l42j: 



■-^12,12 — 



+ 00 



dzdz'(pi{z)(i,2{z)\z - z'|0i(z')02(^') (C2) 



and we recognize the prefactor of (jCl|) to be exactly the 
long wavelength limit of the Coulomb interaction matrix 

element V^'^q ^ : 



~ 2Seeo 

dzdz'cj)xiz)Mz)\z - z'\cj)^,iz')(l)x'iz') (C3) 



Note that we have V^^^^ = V^^q^ = V^^l^ = V^^l\ 
so we recover only the matrix elements that describe the 
interaction of electrons between subbands 1 and 2, as 
should be expected. 

Next, using the fermionic commutation rules, we can 
easily rearrange the binary products of i?-operators into 
four c-operator products: 



l]5|lqS21-q= 4k+q4k'-qClk'Clk (C4) 

q q,k,k' 

X1^21qS21-q= Y cIk+qClk'-qC2k'C2k (C5) 
q qk.k' 

Xl-^21q^21q= ^ 4k+qClk'-qC2k'Clk 



q,k,k' 



H4k+qC2k+q (C6) 
q.k 
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q q,k,k' 

+ cJk+qClk+q (C7) 

q,k 



Note that the sum of the pair terms in (jC6P and ()C7P 
simply acts as the identity operator in the two subband 
subspace and therefore can be ignored. The remaining 
four operator terms can be regrouped in order to provide 
the fong wavelength expansion of the Coulomb potential 
V for our problem: 



Ak+q'^/ik' 



Cu'k'CA'k 



(C8) 



q,k,k' 



The symbol [A/i, fJ.'X'] means that the sum runs only 
along the four Coulomb matrix element mentioned above. 

The final result (jCSp can be interpreted as a dipole- 
dipole interaction, where the oscillating intersubband 
dipole moments interact with each other through their 
local field (sl. 



uj\ - uPl d^A 
2VA dw2 



(D6) 



To transform this equation, we use the following rela- 
tions: 



(D7) 
(D8) 
(D9) 



2 2 2-2 



The first is an immediate corollary from (jD2[) . while 
the second and the third are the Newton formulas for 
equation (jDl[) . We then have: 



2(a;+ -^-f = 2{luI + - 2w+a;_) 



„ _ . dVA 
- 2uJcU}a ) = -TT- 



2K^ 



To obtain the last line, we have derived equation (jD3P 
with respect to w^. This equation clearly leads to the 
result: 



Appendix D: Properties of the polariton dispersion 

The eigenvalue equation for the polariton problem is 
(equation ([ST])): 



= <^ UJr — l^n 



(Dll) 



Furthermore, we can express the minimal splitting: 



/ 2 ~2\/ 2 2\ ro rw 2 2 



(Dl) 



For simplicity we have dropped the wavevector index 
q. The real solution of equation (IDip are obtained from: 



min(w+-^_)2=iii -c^i = (D12) 



This is the result stated in equations (|92|) and (|93|) . 



A =(^2 ^5^)2 _ 4^2-2 

From these equations it is easy to deduce: 



d^ 



4a;+ V dw? / 



The minimal splitting is obtained from: 



= 



(D2) 



(D3) 



(D4) 



(D5) 



Multiplying this equation by uj^ — and using (jD4p 
we obtain: 



Appendix E: Semiclassical plasma Hamiltonian 

In this Appendix we consider the plasma Hamiltonian 
in the semi-classical approximation discussed in the end 
of section IIV CI It will be shown that this Hamiltonian 
leads to the semi-classical dielectric constant (|114p . 

In the semi-classical approximation, it is assumed that 
the intersubband polarization is constant along the het- 
erostructure slab. We can then approximate the current- 
current correlation function by the first order in the ex- 
pansion (jlOOp described in section HVBI 



-i-oo 
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£,a{z)dz / ^I3{z)dz = 
^QW Jo JQ 
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This is equivalent to use a suitably averaged micro- 
scopic response in the quantum well slab. With the use 
of (jEl|l , we obtain the following expressions for the coef- 
ficients that enter the plasma Hamiltonian: 



LcS - Lqw/ fa 

""-Pq ~ Tot 

SSom DTjqw 

CaJ — 1 



(E2) 
(E3) 
(E4) 



We consider first the matter part of the plasma Hamil- 
tonian, which is rewritten as: 



Let A[i jv] be the determinant of iV-th order. Devel- 
oping this determinant along the columns we obtain the 
following recursive expression: 



^[l.-.W] = (^^ - Wi)A[2...Ar] 

N 

—2^ ^ wia;iOiirijiA[2...Ar]/[i] 

i 

-n!2"+^ ^ ujiui.i^...u!i^flu^i}i^i2...ili^^i X 



A 



[2...N]/[n,i2...lr,] 



hujpaUJpp 



baq+P«-q)(p]3-q+P/3q) 



(E5) 



Since the plasmon coupling coefficients are indepen- 
dent from the the index q the latter is dropped in the 
equations, and (IE5I) is expressed in a more handy form: 



= ^ hUapipa + ^ hnapipi + Pa){pl + Pfj) (E6) 



-{N - 1)\2'''UIUJ2-^N^12^23-^NI (ElO) 

with A[2...jv]/[ii,i2 - in] being the determinant of order 
N — 1 — 71, excluding the transitions ii, 12. ..in- The nota- 
tion zi, i2---in means that the sum does not contain repet- 
itive indexes. Using the definition (jE7P the determinant 
is easily rewritten as: 



N 



^[l...N] = {^'^ - Wi)A[2...jv] - ^ '^Pl^^Pi A[2... 



JV]/H 



E-2 -2 -2 A 

^Pl^Pii---^P',^^[2...N]/[n,i2...ir^] 



(E7) 



Let us consider N coupled plasmons; then the Hopfield 
determinant of the resulting Hamiltonian is: 



-{N -l)\Ldl,Ldl^...Ldl^ (Ell) 

Having expressed the A^-th order determinant as a 
function of lower order ones, we make the following re- 
cursive assumption for n < N: 



UJ — UJi 





-ni2 


-ni2 ■■ 


— fllN 


— QlN 





LU + OJ1 


^12^ 


^12 ■■ 


illN 


fllN 






UJ ~ UJ2 





— ^2N 


— ^2N 


fll2 


^12 





UJ + UJ2 ■■ 


^2N 


^2N 



— fllN ^^IN —^2N 
QlN ^IN ^2N 



— ^2N 
^2N 



LO — UJn 

UJ + ujn 
. ^^^) 

To simplify this determinant, we first add every impair 
row to the row above, then we substract every impair 
column from the column on the left. This leads to the 
following simplified determinant: 



UJ — uji 2uj\ 

cj + tJi f2i2 

UJ — UJ2 2tJ2 

ril2 UJ ^UJ2 









^2N 






^2N 



UJ — ujn 2a; at 
UJ + ujj\j 
(E9) 



n n 

A[i...„i=n(-'--')(i-E 



— Ujf 



(E12) 



i=l i=l 

Let us denote Xi — (Iip,j/(cj^ — ujf), then our recursive 
hypothesis implies: 



N 

A[i...7v] - n(a;2 -wf) X 

i=\ 

N N 

(l-a;i)(l-^Xj) -^xia;,(l- ^ 

i=2 1=2 k= ' ' 

N 



N 



k=2,k^i 
N 

k—2.k^ii .12 ■ ..in 

{N -iy.xiX2...XN (E13) 

We can easily show that all the terms in this of sum, 
that imply high orders of the Xi cancel two by two, except 
the linear terms, which completes our recursive demon- 
stration. The final result can be cast in the form: 
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A[i...Ar] 



£qw,ci(w) 



AT 

n( 



(E14) 



where eQw,ci('^) is the semi-classical multi-band di- 
electric constant, provided by equation (|114p . Therefore 
the coupled plasmonic eigenmodes of the Hamiltonian 
(|Et are provided by the zeroes of the dielectric function 
eQW,ci(w). 

If we now add the light-matter interaction with a single 
photonic mode, similar reasoning leads to the Hopfield 
determinant: 



A(c.) 



eQw,ci('^) 



N 



2 2 



£oft,cl(t^) 



(E15) 



with £cff,ci(w) the classical effective medium constant 
defined by: 



l-/u 



eoff,ci(w) eQw,ci(w) 



(E16) 



Here fw = Lqy^ / Lcav is the filling factor of the quan- 
tum well (heterostructure) in the cavity, which is inde- 
pendent from the intersubband transition in this approx- 
imation. 

From a mathematical point of view, the analytical di- 
agonalisation of the problem was obtained thanks to the 
relation Ca^p = 1- This would not be possible in the gen- 
eral case, and yet it is clear that the Hopfield determinant 
of the Hamiltonian (1551) is a polynomial of 2 x [N + l)-th 
degree. Therefore we can still use the definition (I107P to 
compute numerically the exact quantum effective dielec- 
tric constant of the problem. 
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